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In classical computational chemistry, the coupled-cluster ansatz is one of the most commonly used 
ah initio methods, which is critically limited by its non-unitary nature. The unitary modihcation as 
an ideal solution to the problem is, however, extremely inefficient in classical conventional computa¬ 
tion. Here, we provide the hrst experimental evidence that indeed the unitary version of the coupled 
cluster ansatz can be reliably performed in physical quantum system, a trapped ion system. We 
perform a simulation on the electronic structure of a molecular ion (HeH^), where the ground-state 
energy surface curve is probed, energies of excited-states are studied and the bond-dissociation is 
simulated non-perturbatively. Our simulation takes advantages from quantum computation to over¬ 
come the intrinsic limitations in classical computation and our experimental results indicate that 
the method is promising for preparing molecular ground-states for quantum simulation. 

PACS numbers: 03.67.Ac, 31.15.Dv, 37.10.Ty, 42.50.Dv 


The central problem in quantum chemistry and molec¬ 
ular physics is to determine the electronic structure and 
the ground-state energy of atoms and molecules by solv¬ 
ing the quantum many-body equations, which is gener¬ 
ally intractable due to the exponential scaling to the size 
of the system. Quantum simulation [IHS] can provide 
the solution for such “exponential catastrophe” prob¬ 
lem. The key ingredient of quantum molecular simulation 
consists of (i) ground (excited) -state preparation and 
{a) energy estimation of the corresponding state mill. 
Recently, the assessed costs for the energy estimation 
for a well-prepared ground-state in quantum computa¬ 
tion have been immensely reduced CHII], indicating that 
chemistry simulation can be one of the main applications 
of a quantum computer in near future. However, it is still 
remaining major obstacle to efficiently and reliably find 
the molecular ground state, which belongs to the class of 
extremely hard problems called Quantum Merlin Arthur, 
the quantum analog of NP-hard problem [mils]. Re¬ 
cently various theoretical schemes for the ground-state 
problem have been proposed and proof-of principle ex¬ 
perimental demonstrations have been performed includ¬ 
ing the adiabatic pAHlG] and algorithmic preparations 

[I7H20]. 

For the ground-state problem, the developments of 
conventional quantum chemistry can be adopted to quan¬ 
tum computation. In computational chemistry, it has 
been the main focus to circumvent the problem by ap¬ 
proximating the many-body Schrodinger equation and a 
series of theoretical and numerical methods have been 
developed. The coupled-cluster method is one of the 
most prominent ab initio methods for finding a molec¬ 
ular ground state and it is considered to be the cur¬ 
rent gold standard [HHH]. However, the coupled-cluster 
ansatz is built with non-unitary operation, which leads 
to drawbacks such as lacking a variational bound on 


the ground-state energy [221126] . The unitary version 
of the coupled-cluster methods would perfectly resolve 
the problem, whereas it is classically inefficient without 
proper truncation of the infinite series expansion. It has 
been a long-standing challenge to build an efficient com¬ 
putational scheme for the unitary coupled-cluster (UCC) 
ansatz. The authors of Refs. 0 HZ] pointed that the 
UCC ansatz can be efficiently implemented in a quantum 
computer. In other words, the quantum implementation 
of the UCC method can outperform the classical compu¬ 
tation for the problem of finding molecular ground-state. 

Recently a variational method of approximating molec¬ 
ular ground states has been experimentally tried in a pho¬ 
tonic system HZ). Due to challenges in the system, how¬ 
ever, the variational ansatz employed in the experimental 
demonstration was not the UCC ansatz, but called ’de¬ 
vice ansatz’ which is device specific method, therefore is 
not scalable. As discussed in Ref. [6], the UCC ansatz 
provides a generic and scalable scheme for generating a 
parameterized state for the variational method and can 
be implemented efficiently with quantum devices includ¬ 
ing trapped ions. 

Here we report the first experimental realization of 
the UCC ansatz with a minimal basis, based on quan¬ 
tum simulation in a multi-level of a trapped ^^^Yb+ ion. 
We simulate the electronic structure of a molecular ion 
(HeH+) [27l EHl and reliably find the molecular ground- 
state as well as the corresponding energy by the UCC 
ansatz and the variational method, which can be consid¬ 
ered as an alternative method for the energy estimation 
[291130]. Moreover, we apply the quantum UCC method 
to compute excited states and chemical-bond softening 
non-perturbatively. 

The coupled-cluster approach is based on the varia¬ 
tional method with the trial state of the exponential 
ansatz in the form of \G). Here \G) is a reference 
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FIG. 1. (a) The conceptual procedure and (b-f) the experi¬ 

mental realization of the UCC simulation to find the ground 
state of target molecule, (b) The Hartree-Fock basis states 
of the target molecule, HeH^. (c) The mapping of the basis 
states on the energy levels of ^^^Yb^ including cluster ampli¬ 
tudes and t 2 , which are controlled by the duration of 

microwave pulses, (d) The microwave pulse sequence for the 
preparation of the UCC ansatz. The effective time evolution 
operator is expanded by the Suzuki-Trotter scheme 

(see SM. C: Quantum Implementation of the UCC Ansatz). 
(e) The measurement of energy (H) of the UCC ansatz given 
cluster amplitudes, (f) A classical minimum search algorithm 
(Nelder-Mead) is applied to determine the cluster amplitudes 
for the next UCC ansatz and the final ground state. 

state, such as the Hartree-Fock ground state and the 
cluster operator T is constructed by a sum of n-electron 
excitation operator with the transition amplitudes as 
variational parameters (see Eq. (SI) in Supplemental Ma¬ 
terial (SM)). However, the operator T is not necessarily 
Hermitian and energy estimation by the ansatz is not 
guaranteed to bounded by variational theorem [25l [26] . 

The UCC scheme based on the form of the following 
ansatz 

I'i/’ucc) = |G) (1) 

which apparently provides a solution of the non- 
Hermitianity problem in the coupled-cluster theory. 


However, classical implementation of the UCC have in¬ 
trinsic limitation, e.g., infinite series of the expansion 
[26] (see also SM. B: Classical Implementation of UCC 
Ansatz). As a result, all classical applications of UCC 
involve some type of truncations with potentially uncon¬ 
trollable errors. On the other hand, the unitary operator, 
U = , can be considered as a time-evolution opera¬ 
tor, Te., U = driven by an effective Hamiltonian 

= i(T — T”*") with a dimensionless time interval set 
to be 1. Since the time-evolution is efficiently simulated 
in a quantum system m, the quantum implementation 
of the UCC ansatz can reduce the computational cost 
much less than the classical requirement. 

The whole procedure of finding the ground state of 
a target molecule is shown in Fig. [^a), which is also 
discussed in Refs. 0117 ]. After efficiently preparing a 
trial state with UCC ansatz in a quantum system that 
maps the classical basis set of the target molecule, we 
measure the average energy of the state. The prepara¬ 
tion of the UCC ansatz and the energy measurement are 
performed in the quantum system. Based on a classical 
feedback algorithm, we adjust the parameters, i.e., the 
cluster amplitudes of the UCC ansatz. We repeat the 
quantum process of preparation and measurement until 
we find the variational minimum of the target Hamilto¬ 
nian. 

Target Molecule (see also SM. C: The Hamiltonian 
of HeH+). We choose the helium hydride cation (HeH+) 
m [28] for the computation of the energy curve. In 
the second quantization representation with the minimal 
STO-3G basis-set [24] from the Is orbitals of Hydrogen 
and Helium, the Hamiltonian of HeH+ is described by 

H{R) = E hpq {R)alaq + I E ^pqrs i^lVjOj^rpOj^qCLr^s (2) 

pq pqrs 

where R is the nuclear separation between hydrogen and 
helium, hpq{R) and hpqrs{R) are related to one electron 
and two electron transitions, respectively and the index 
p^q^r^s stands for the four possible states in our Hilbert 
space. The terms of hpq{R) and hpqrs{R) are computed 
numerically by the Hartree-Fock method with the scaling 
of 0{M^)^ where M is the number of molecular orbits. 
The creation and annihilation operators in the Hamil¬ 
tonian (§ are mapped to spin Pauli operators by per¬ 
forming the Jordan-Wigner transformation and pairs of 
Pauli operators are mapped to four-level systems. After 
the mapping, the Hartree-Fock basis for HeH+ consists of 
the following set of four states, {|G), \Eii) , \Ei 2 ) , IE 2 )} 
as shown in Fig. [^b). 

Mapping of HeH+ on ^^^Yb^ ion. The electron 
excitation operators, which excite the electrons out of the 
Hartree-Fock ground state, up to two electron excitations 
are given by, 

Ti = + ti2a2^ai^ , T 2 = t2a\^al2^ai^aii (3) 
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Note that all the terms are spin preserving, and the 
^ 115 ^ 12,^2 are in general complex numbers to be deter¬ 
mined by an optimization process. After the same map¬ 
ping process of the Hamiltonian ([^, the effective Hamil¬ 
tonian i^eff = i(T — T’f') for the cluster operators is writ¬ 
ten as 


= ^tn(|En)(G| + |E2)(Ei2|) 

+ iti2 (1^12) (G| +1^2) (EllI) (4) 

+ it 2 IE 2 ) (G| + h.c.. 

We realize the effective Hamiltonian in a quantum 
system of multiple energy levels in trapped ^^^Yb^ ion. 
As shown in Fig. Re), four energy levels in the ground- 
state manifold of ^ 1/2 of the ^^^Yb+ are employed [32l 
[ 33 ] to map the basis state as |E = 0, = 0) = |G) and 

\F = l,mF = -1,1,0) = {|Eii), IE 12 ), IE 2 )}, which are 
separated by cjhf — , ^hf + and cjhf, where the 

hyper-fine splitting of cjhf = (27r) 12.642821GHz, Zee- 
man splitting of ujz = (27r) 13.586MHz with the static 
magnetic field of E = 9.694G. 

Preparation of the UCC Ansatz (see also SM. D: 
Quantum Implementation of the UCC Ansatz). The uni¬ 
tary operator U = is implemented as a time evolu¬ 

tion of the system with the effective Hamiltonian Eeff as 
shown in Fig. Rd). The initialization of the state to |G) 
is performed by the standard optical pumping technique. 
The transitions {|G) jEn), |G) IE 12 ), |G) IE 2 )} 

are implemented by applying resonant microwaves. 
The other transitions {|Eii) IE 2 ), IE 12 ) ^ IE 2 )} are 
achieved by applying composite pulse sequences shown 
in the insets of Fig. Bd). Consequently, the experimen¬ 
tal implementation of the unitary operator U is achieved 
by the sequence depicted in Fig.Rd), which results from 
the second-order Suzuki-Trotter expansion. In the exper¬ 
iment, the UCC amplitudes tn, ti 2 , and ^2 are controlled 
by the durations of the corresponding microwave transi¬ 
tions. We note that for HeH+, the amplitudes near the 
molecular ground state are much smaller than 1 and the 
errors from small Trotter expansions {N = 2 in our ex¬ 
periment) are negligible. 

Energy Measurement of the UCC Ansatz. We 

can obtain the energy {H) = Y,pq +Ep 9 rs i^pgra), 

where Hpq = Upga^Uq and Hpqrs = by 

term-by-term measurements and addition of all of them 
in the target Hamiltonian (R. For our case of HeH+ in 
the minimal basis, it requires 24 times of measurements 
(SM. C. The Hamiltonian of HeH+), which necessarily 
needs the information of all the components in the den¬ 
sity matrix of the UCC ansatz (see also E. Measurement 
of the Energy for HeH+). Note that as the system size 
increases, we do not need the full knowledge of the den¬ 
sity matrix of the state for the energy measurement, since 
the number of terms in the Hamiltonian <§ scales poly- 
nomially IS [27]. Since we need the full knowledge of 
the density matrix for our small scale simulation, we re- 



FIG. 2. The search process of the minimum energy at 
R = 1.7a.u. assisted by the classical Nelder-Mead algorithm 
with UCC ansatz. The measured energy (H) (dots) and the 
fidelity of the prepared state (bars) to the ideal ground state 
depending on the number of iterations (a) with full six pa¬ 
rameters and (b) with two parameters. For both cases, the 
algorithm converges to the ground-state energy obtained by 
the exact diagonalization with decent fidelity of the state. 
Red dots show the successful steps that contribute to the 
convergence, (c) The side view and (d) the bottom view of 
the searching process with two parameters for the successful 
steps. The atomic unit (a.u.) is used for all the figures. 

construct the full density matrix by the standard 

quantum state tomography, which requires 15 times of 
measurement, and obtain the energy by Tr For 

the relevant components of the density matrix, we repeat 
the standard measurements up to 1000 times, which give 
3.2% projection uncertainty of standard deviation. 

Classical Minimization Algorithm The prepara¬ 
tion and the measurement of an UCC ansatz are per¬ 
formed in quantum system and the minimization process 
is supported by classical algorithm. The measured value 
of {H) for the prepared UCC ansatz is taken as an in¬ 
put for a classical optimization algorithm, which com¬ 
pares it to the previous values and suggests a new set of 
{Ui, U 2 , ^ 2 } so that the same procedure is repeated until 
the resulting (H) converges to some value. As a result, 
we obtain an optimized state with minimal energy for 
approximating the ground state of HeH+ in the form of 
an UCC ansatz in Eq. 0 . In our realization, we use a 
popular Nelder-Mead minimum search algorithm [34]. 

Fig. 15 shows an instance of the energy optimization 
process, when the nuclei separation of HeH+ is fixed to 
be R = 1.7 a.u.. Note that throughout the paper, the 
atomic unit (a.u.) is used. The algorithm is capable of 
finding the minimum energy and state in around hun¬ 
dred iterations with the full six-parameter simulations as 
shown in Eig. Ra). About twice less iterations shown 
in Eig. Rb) can be achieved for an ansatz simplified to 
contain two parameters (see SM. E: Reduction of Pa- 
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Nuclei distance R 



FIG. 3. The ground state energy of HeH^ depending on the 
inter-nucleus distance R. The error-bars of the experimental 
data mainly come from the quantum projection noise of 1000 
repetitions for each term of the Hamiltonian §. 

rameters). Since both cases provide equivalent results, 
we focus on the two-parameter ansatz in the following 
discussion. Figs.|^c)(d) show the typical search of mini¬ 
mum energy by the classical Nelder-Mead algorithm with 
two parameters. 

Fig-i shows the energy curve of the ground state of 
HeH+ depending on the nuclear distance R, where each 
point is obtained by the procedure of Fig. The exper¬ 
imental data are in agreement with the energy (orange 
line) calculated by the exact diagonalization of the full 
matrix of Hamiltonian within the error bars. From 
the energy curve, the equilibrium distance between the 
nuclei is located at i? = l.TSa.u. with the corresponding 
energy of E = —2.86 ± 0.05 a.u.. 

Furthermore, the same procedure can be used to study 
the non-perturbative behaviors of the HeH+ molecu¬ 
lar ion under strong electric field with the new target 
Hamiltonian including the effect of the electric field as 
E • (ri + r 2 ) — E • (2RHe + Rh) (see SM. F: The Electric 
Field Effect on HeH+ ). Eig.j^a) shows the phenomenon 
of chemical-bond softening of HeH+ (at R = 1.7 a.u.) as 
the strength of the electric field increases, which even¬ 
tually leads to a dissociation of the molecular ions [35] . 
We compare our non-perturbative results with those ob¬ 
tained through the first-order and second-order pertur¬ 
bation theories shown in Eig. [^b). 

Einally, we also study the excited states of H by chang¬ 
ing the target Hamiltonian to (i7 — A)^, where A is a 
parameter close to the energy of an excited state, which 
turns the excited state of H into the ground state of 
(i7 — A)^ (see SM. G:The Computation of Excited States 
Energy). In the experiment, we uniformly scan the val¬ 
ues of A and apply the same UCC procedure to find the 
minimum energy in a given A. As shown in Eig. gc), we 
observe that the required precision for the computation 
of excited states should be much higher than the sepa¬ 


FIG. 4. Applications of the UGG simulation, (a) The ground 
state energy of HeH^ subject to a static electric field along 
the nuclei axis for different strengths, (b) The comparison 
between the UGG quantum simulation and the perturbation 
theory at given inter-nucleus distance R = 1.7a.u.. (c) The 
search process of the energies of excited states of H by find¬ 
ing the ground-state energy of the Hamiltonian {H — A)^ by 
scanning the values of A. When A is the energy of an excited 
state, the experimental minimum value of ((if — A)^) tends to 
be zero. We also can calculate the excited-state energy from 
other non-zero values of ((if — A)^). If A is on the left (right) 
side of the excited energy, the positive (negative) solution of 
Emeas = ((if — A)^) provides the excited energy. 

ration of the energies. In the current limited system, we 
obtain the energy of the highest excitation that has rel¬ 
atively large energy gap to other states but the rest of 
them are not well resolved. 

Our current realization is capable of simulating any 
molecule up to four electronic levels with a single ion. In 
general, a molecule of N electrons system in M molecu¬ 
lar orbitals (M > N) can be implemented with M qubits 
system or M/2 qudits, four-level systems shown in our 
realization, through the Jordan-Wigner transformation 
and four-level mapping. Eor the UCC implementation 
with M qubit system, it requires the simulation of time- 
evolution of M-body interaction, which is equivalent to 
the nonlocal product of M Pauli operators. The sim¬ 
ulation of such M-body interaction, which is the most 
challenging operation in the UCC protocol, can be per¬ 
formed by applying 2M times of CNOT-gate or 2 times of 
the multi-particle M0lmer-S0rensen gates [HIMlEl]. The 
measurement of M-qubit Hamiltonian with the O(M^) 
terms has been already well established in the trapped ion 
system. Eor the M/2-qudit Hamiltonian, we can simply 
use the same measurement scheme used in our experi¬ 
mental demonstration. The UCC scheme for the trapped 
ions can be applied to other physical platforms [381440] . 

We emphasize the computational complexity of the 
quantum implementation of the UCC method scales 
polynomially with the number of orbitals M. Including 
the maximum excitation up to f, each cluster operator 
contains k creation operators and k annihilation opera- 
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tors. For a total of M orbital modes, therefore, we have a 
total of terms. After the Jordan-Wigner trans¬ 

formation, the fermionic operators are mapped into spin 
operators, which requires 0{M) operations. The total 
number of scaling as the number of molecular orbits M 
is Moreover, the time evolutions and the 

measurements in our UCC implementation allow parallel 
computation ISIIT], which boosts the performance. Our 
experimental realization of UCC method opens a new 
dimension of quantum simulation and offers a solution 
for the classical coupled-cluster methods. We note that 
some of other current developments and understandings 
in the coupled-cluster schemes could be adapted in quan¬ 
tum UCC scheme. Moreover, our UCC scheme could be 
applied to other large eigenvalue problems in network 
search algorithm and condensed matter physics. 

Notes This paper is a revised version of 
arXiv: 1506.00443 by the same authors. While the 
work is under consideration for publication, the authors 
noted that a paper reporting Scalable Quantum Simula¬ 
tion of Molecular Energies in Phys. Rev. X 6, 031007 
(2016). 
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